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Abstract 

While modern optics is largely a physics of harmonic oscillators 
and two-by-two matrices, it is possible to learn about some hidden 
properties of the two- by- two matrix from optical systems. Since two- 
by-two matrices can be divided into three conjugate classes depending 
on their traces, optical systems force us to establish continuity from 
one class to another. It is noted that those three classes are equivalent 
to three different branches of Wigner's little groups dictating the in- 
ternal space-time symmetries massive, massless, and imaginary-mass 
particles. It is shown that the periodic systems in optics can also be 
described by have the same class-based matrix algebra. The optical 
system allow us to make continuous, but not analytic, transitions from 
massiv to massless, and massless to imaginary-mass cases. 

1 Introduction 

Two-by-two matrices with real elements have three independent parameters 
if their determinants are constrained to be one. They constitute building 
blocks for many branches of physics, including beam-transfer matrices in 
optics p] and Wigner's little groups for internal space-time symmetry of 
particles [2 EH 13]. 

For the two-by-two matrix, we are accustomed to solve a quadratic equa- 
tion to get the eigenvalues and construct a rotation matrix to get the eigen- 
values. This procedure does not always lead to correct answers, because 
squeeze matrices should also be considered [S]. We are quite familiar with 
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rotations, but the concept of squeeze started getting our attention only after 
squeezed states of light appeared in the physics literature [5] . 

In particle physics, Lorentz boosts are squeeze transformations, and this 
aspect was addressed by Paul A. M. Dirac in his 1949 and 1963 papers [8|. 
It is possible to apply this concept to high-speed hadrons which are bound 
state of the quarks which are thought to be more fundamental particles [3J. 
Thus, high-energy hadronic physics and modern optics share the same math- 
ematical base, and it is profitable to trade physics between these two branches 
of physics using the common mathematics language. 

In this report, we note first that optical beam transfer matrix, often called 
the ABCD matrix, has three independent parameters and its determinant is 
one. We study then how this matrix can be decomposed into one-parameter 
matrices. The Bargmann decomposition and Iwasawa decomposition are 
already familiar to us, and are used often in the literature |9] [10] . 

We show that, in addition, there is a decomposition based on the concept 
of conjugate classes [TTJ. There are three conjugate classes. The first class 
consists of those matrices with their traces smaller than two, the second class 
consists of those with the traces equal to two, and the third consisting of 
those matrices with their traces greater than two. It is remarkable that this 
purely mathematical theorem corresponds to Wigner's construction of his 
little groups, which dictate the internal space-time symmetries of massive, 
massless, and imaginary-mass particles respectively 

In Sec. HI we introduce the Wigner decomposition based on the conjugate 
classes of the two-by-two matrices. It is then shown that the Wigner decom- 
position can be translated into the Bargmann and Iwasawa decompositions. 
In Sec. |3j we discuss how we can formulate those three decompositions while 
studying optical multilayer systems. It is shown that the Wigner decompo- 
sition is needed for repeated application of the ABCD matrix for periodic 
system. In Sec. HJ we study how those three decompositions can serve useful 
purposes in studying Wigner's little groups and thus the internal space-time 
symmetries of elementary particles. 



2 Decompositions of the ABCD Matrix 



The two-by-two matrix with real elements and unit determinant can be writ- 
ten as 




(1) 
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with AD — BC = 1 has three independent elements. This matrix is commonly 
used as the beam transfer matrices in optics [T]. The complete set of these 
matrices is like the group SU(1, 1) which serves as the fundamental language 
for squeezed states of light [H [12] 

All the matrices in this set can be divided into three classes depending 
on their traces [11J. Without changing its trace, we can bring every matrix 
to an equi-diagonal form by a rotation |15j . 

We shall use the notation [ABC A] for the equi-diagonal ABCD matrix. 
This equi-diagonal matrix now has two-independent parameters. If its trace 
is smaller than two, the matrix can be written as 

\ABCA] — ( COS( ^ - er?sin <^ ( 2 ) 
\ e~ v sin <p cos 4> ) ' 

If the trace is greater than two, it can take the form 

\ABCA}=( C ° shx -e v sinh X \ ^ 
\—e~ v smhx coshx /' 

If the trace is equal to two, it can be brought to the form 

[ABCA] = (l -f). (4) 

. We choose to use the notation W(r) collectively for the following three 
matrices. 

/cos0 — sin0\ / coshx — sinh^A / 1 ~l\ /r\ 
V sin (f) cos <fi ) ' \ — sinh x cosh \ ) ' V 1 / ' 

The parameter r could be 8, x, or 7. If we define the matrix B{rf) as 

Then every equi-diagonal ABC A matrix can be written as 

[ABC A] = B(r])W(T)B(-r]). (7) 

Indeed, the matrix W(t) constitutes a set of three matrices which play an 
important role in the theory of two-by-two matrices. 
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It is important to note that these three matrices constitute the basic ele- 
ments for Wigner's little group which dictates the internal space-time symme- 
tries of elementary particles [2],|3]. We thus choose to call W(t) the "Wigner 
matrix" [13] and call the expression of Eq.([7|) "Wigner decomposition." The 
Wigner decomposition leads to 

[ABCAf = B(rj)[W(r)} N B(-r]) = B(t])W(Nt)B(-t]), (8) 

convenient in dealing with repeated applications of the ABC A matrix in 
periodic systems. 

It is known that the ABC A matrix can be written as the product of three 
one-parameter matrices 

cos(#/2) sin(0/2) \ / cosh A -sinh A \ / cos(#/2) sin(0/2) \ 
-sin(0/2) cos(0/2)J U sinh A cosh A ) \-wi(9/2) cos(0/2) J ' 

(9) 

This form is called the Bargmann decomposition j£l [H] . This expression can 
be compressed to 

r cosh A cos 9 — sinh A — cosh A sin 9 \ , , 

\ — sinh A + cosh A sin 9 cosh A cos 9 J ' 

with two independent parameters. We shall call this expression the "Bargmann 
matrix." 

If sin# = tanhA, the ABC A matrix takes the triangular form given in 
Eq.flSJ). Then this special case of the Bargmann decomposition is called the 
"Iwasawa" decomposition. This form is also known to correspond to gauge 
transformations of massless particles [31 Q3] . 

We shall see in this report that both Wigner and Bargmann decompo- 
sitions play essential roles in optics and space-time symmetries. The ques- 
tion then is whether one transformation can be translated into the other. 
If cosh A cos 9 is smaller than one, the off-diagonal elements have opposite 
signs. The Bargmann matrix of Eq. fllOl) will be translated into Eq.(J2]) with 



cos d> = cosh A cos 9, eP — \ — — - — , ^ (11) 
Y V sin^ + tanhA v ; 

If the diagonal element is greater than one, the matrix will be translated into 
Eq.fl3D with 



/tanhA — sm 9 , , 

cosh x = cosh A cos 9 , e 1 = \ — (12) 

V tanh A + sm 9 
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Figure 1: Transitions from sin to sinh, and from cos to cosh. They are 
continuous transitions. Their first derivatives are also continuous, but the 
second derivatives are not. Thus, they are not analytic continuations. 



These variable transformations are simple enough, but the real issue is what 
happens when the diagonal element makes a transition from less-than-one to 
greater-than-one. If it becomes one, the lower left element of the Bargmann 
matrix becomes zero, and the matrix becomes triangular. 

In order to study this transition, we introduce a small number 



and see what happens when it changes its sign. When e is positive, the 
Bargmann matrix can be written as 



e = cosh A sin 6 — sinh A, 



(13) 



( 



1 — e sinh 77 cosh r] 
e cosh i] 



—2 sinh 7] 
1 — e sinh 77 cosh r] 



) 



(14) 



If we let 




(15) 



the matrix becomes 




(16) 



which can be decomposed into 




(17) 



If e is negative, we should define a and (3 as 



and the matrix becomes 



(l + a 2 /2 
V -a//3 



—a/3 
1 + a 2 /2 



) 



(19) 




1 + a 2 /2 



) 



( 



VP o 





(20) 



It is clear now that the transition from the rotation-like Wigner (with 
cos0 as the diagonal element) to the squeeze-like matrix (with coshx as 
the diagonal element) is like the transitions from sin a to sinha, and from 
cos a to cosh a, as shown in Fig. [TJ This is a continuous transition, but not 
analytic. The second derivative is not continuous. 

3 Periodic System in Optics 

Let us consider an optical beam going through multiple layers consisting of 
two different refractive indexes. This problem has been extensively discussed 
in the literature. The two-by-two matrix formulation of this problem is given 
in the Appendix. 

When the beam goes through the first medium, we can use the matrix 



For the second medium, we use 02 instead of 0i. 

If the beam in the first medium hits the second medium, it is partially 
transmitted and partially reflected. According to the Appendix, the bound- 
ary matrix takes the form 



This form is also given in Eq.([6]) in connection with the Wigner decompo- 
sition. When the beam hits the first medium from the second medium, the 
boundary matrix is B(—rj). 

Let us consider the cycle which starts from the half-way in the second 
medium and ends at the second medium as illustrated in Fig [2J Then the 
beam transfer matrix becomes 




(21) 




(22) 



[ABC A] = R (<f> 2 ) [B(q)R (200 B(-r,)) R (0 2 ) . 



(23) 
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Figure 2: Optical layers. There are phase-shift matrices for their respective 
layers. There is a boundary matrix for the transition from the first to second 
medium, and its inverse applies to the transition from the second to first 
medium. The cycle starts from the middle of the second layer. 

The quantity inside the square bracket is takes a form the Wigner decom- 
position, but it is sandwiched between two rotation matrix. Our problem is 
to write the entire matrix chain as a Wigner decomposition, convenient for 
the periodic system. For this purpose, we write the Wigner decomposition 
in the middle as a Bargmann decomposition 

R{6)S{-2\)R{6), (24) 

where 9 and A are determined by 



cosh A = (cosh rj) \l 1 — cos 2 0i tanh 2 r\ 



a COS (pi 

cos# = , . (25) 



(cosh rj) \J 1 — cos 2 0i tanh 2 r) 

We can now write [ABC A] as 

R{<p 2 )R{0)S{-2X)R(d)R{<p 2 ). (26) 

Since R(8)R (02) = R(9 + 02/2) , the [ASCA] matrix becomes a Bargmann 
decomposition of the form 

[ABC A] = R(6*)S(-2\)R(e*), (27) 
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with 

9* =9+ Ifa. (28) 

After matrix multiplications, the [A.BCA] matrix can be written as 

/ cosh A cos 9* — sinh A — cosh A sin 9* \ , . 

\ — sinh A + cosh A sin 9* cosh A cos 9* J 

For the multilayer system, we need a repeated application of this from. 
For this purpose, we have to convert this into the Wigner decomposition: 

B{rf)W{T*)B{-rf). (30) 

When cosh A cos 9* is smaller than one, 

\xr(«*\ f cos< ^* -sin0*\ , . 

W(r)= Uin0* cos0* J' (31) 

with 



j* i a n* f *\ sin 9* — tanhA , nn . 

cos = cosh A cos 9 , exp (77 ) = \ . Qt . , t-t (32) 

V sm 9* + tanh A 

If the diagonal element is greater than one, the W{r*) matrix is of the form 

W(r*)=( C ° S \ X \ - Si f **Y (33) 
V — sinh x cosh x J 

with 



/ tanh A — sin 9* 

cosh x = cosh A cos 9 , exp (rj) = ^ ^^—^ (34) 

The situation is similar if the diagonal element is one, and the matrix is 
triangular. We can use this Wigner decomposition to calculate 

[B{t]*)W ( t*)B{-t]*)] n = B{ V *) [WW*)] N B(-ij) (35) 

It is seen that this periodic system contains the Wigner, Bargmann, and 
Iwasawa decompositions in its natural language, and it forces us to bring the 
[ABCA] matrix into its Wigner decomposition for its repeated applications in 
the periodic system. In SecHJ we shall study what this Wigner decomposition 
means in symmetries in particle physics. 
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4 Space-time Symmetries 



In 1939 [2], Eugene Wigner considered subgroups of the Lorentz group whose 
transformations leave the given momentum of a particle invariant. These 
subgroups are called Wigner's little groups. While leaving the momentum 
invariant, the little group transforms internal space-time variables. For in- 
stance, if the particle is at rest, rotations do not change its momentum, but 
they can change the orientations of the particle spin. In this section, we for- 
mulate Wigner's little group using the technique of decompositions discussed 
in Sec. |2]and Sec. |3j 

In Sec. El we studied the decompositions of the [ABCD] into single- 
parameter matrices, using the two-by-two matrices B(rf), R(<f>), and S(x)- 
It is known that these two-by-two matrices correspond to the four-by-four 
Lorentz-transformation matrices applicable to the Minkowskian space-time 
coordinate variables (z,y,z,t) [3j. The two-by-two matrix B(rj) can also 
be written as the the four-by-four matrix for the Lorents boost along the z 
direction: 



B( v ) 



-rill 






Vo 





1 









cosh 7/ 
sinh 7/ 






sinhr/ 
cosh T] J 



(36) 



which performs a Lorentz boost along the z direction. The two-by-two matrix 
R{4>) an be translated into 



R{4>) 



cos(0/2) 
sin(0/2) 



sin 



K0/2) \ , 
cos(0/2) ) 



/ COS( 





V 



sin i 




sin < 

1 
cos < 







1/ 



(37) 



which performs a rotation around the y axis. The S(x) corresponds to 
'cosh(x/2) sinh(x/2)\ 



six) 



sinh(x/2) cosh(%/2) / 



/ coshx 



\ sinh x 




1 








1 





sinh x \ 



cosh x J 



which performs a Lorentz boost along the x direction. 

It is possible now to translate the contents of Sees. |2] and |3] into the 
language of four-by-four Lorentz transformation matrices applicable to the 
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Minkowski space of (z,y,z,t). In this convention, the momentum-energy 
four- vector is (p x ,p y ,p z , E) . If the particle moves along the z direction, this 
four- vector becomes 



0,0,p,^ 2 + m 2 J, (39) 

in the unit system where c = 1, where m is the particle mass. We can obtain 
this four-vector by boosting a particle at rest with the four-momentum 

(0,0, 0,m), (40) 

using the four-by-four boost matrix given in Eq. fl36|) . with 

p_ 

y/p 2 + m 2 



tanh^) = , / (41) 



Now the four-momentum of Eq. (140]) is invariant under the rotation matrix 



R{2<P) 



( cos(20^) sin(2(/>) 0\ 

10 

-sin(20) cos(20) 

V o 001/ 



(42) 



Thus, the matrix 

B{r ] )R{2 ( j ) )B{-ri) (43) 

leaves the four-momentum of Eq. fl39|) invariant. After making this rotation, 
we can bring the momentum to its initial state by boosting it by B(rj). The 
net effect is the momentum-preserving transformation. This set of transfor- 
mations is illustrated in Fig. El and corresponds to the Wigner decomposition. 

If the particle has a space-like momentum, we can start with the four- 
momentum 

(0AP,E), (44) 

where E is smaller than p, which it can be brought to the Lorentz frame 
where the four-vector becomes 

(0,0,p,0). (45) 
The boost matrix takes form of Eq. (l36|) . with 

tanh(r?) = — . (46) 
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M X 



Rotate without 
changing momentum 



Boost 



Momentum 




Figure 3: Illustrations of the Wigner decomposition (left) and the Bargmann 
decomposition. In both cases, the net transformation leaves the four- 
momentum of the particle unchanged. 



The four- momentum of Eq.( H5|) is invariant under the boost 

( cosh(2x) 



V - sinh(2x) 

along the x direction. 

Let us finally consider a massless particle with it four-momentum 



S(-2 X ) 



- sinh(2x) \ 

1 
1 
cosh(2x) / 



(47) 



(0,0,p,p). 



(48) 



/ 1 


27 
V2 7 




1 





(49) 



It is invariant under the rotation around the z axis. In addition, it is invariant 
under the transformation 

-27 27 

1 - 2 7 2 2 7 2 
-2 7 2 l + 2 7 2 / 

This four-by-four matrix has a stormy history [31 0] , but the bottom line is 
that it corresponds to the triangular matrix of Eq.([5]), and the variable 7 
performs gauge transformations. 
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We can obtain this massless case from the massive or imaginary case using 
the limiting procedure spelled out in Sec. [2J This procedure is widely known 
as the group contraction in the literature. 

As for the Bargmann decomposition, let us go to Fig. [3j The particle 
moving along the z direction can be rotated first by R(6). It can then be 
boosted along the negative x axis, and then rotated again by R(9) to the 
original position. Indeed, this is a momentum-preserving transformation. 
The net transformation can be written as 

R{0)S(-2\)R(0). (50) 

This Bargmann decomposition is applicable to all three cases of the momen- 
tum inj. 

In addition, Wigner's little group allows rotations around the momentum 
which does not change it. This extra degree of freedom does not affect the 
description of the symmetries given in this section [IB] 

Concluding Remarks 

In this report, we started with a two-by-two matrix with real elements, but we 
were led to consider three classes of equi-diagonal matrices, with their traces 
less than two, equal to two, and greater than two. From the mathemati- 
cal point of view, this process is the construction of group representations 
according to conjugate classes. 

These conjugate classes correspond to Wigner's little groups for the inter- 
nal space-time symmetries for massive, massless, and imaginary-mass parti- 
cles in the Lorentz-covariant world. Indeed, this aspect is a very "unreason- 
able agreement" between mathematics and physics. 

It was noted that these equi-diagonal matrices as products of three one- 
parameter matrices, resulting in the Wigner, Bargmann, and Iwasawa de- 
compositions. It is noted further that the optical periodic systems, such 
as multilayer optics, can perform these decompositions. Thus, the optical 
periodic system speaks the language of the fundamental symmetries for ele- 
mentary particles in Einstein's Lorentz-covariant world. 
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Appendix 



The iV-layer optics starts with the boundary matrix of the form [T7J [TO] 

n(n\ — cosh(?7/2) sinh(?7/2)\ f . 

{V) Uinh(V2) cosh(V2) J ' { > 

which, as illustrated in Fig. [21 describes the transition from medium 2 to 
medium 1, taking into account both the transmission and reflection of the 
beam. As the beam goes through the medium 1, the beam undergoes the 
phase shift represented by the matrix 

P(<h) = { JL)- (52) 

When the wave hits the surface of the second medium, the corresponding 
matrix is 

B(-V) = ( C0S \ ( ? / 1 " S ^ (r !{? ) , (53) 



— sinh(?]/2) cosh(r//2) / 

which is the inverse of the matrix given in Eq. (15 lj) . Within the second 
medium, we write the phase-shift matrix as 

W 2 ) = ( 6 T e °0" (54) 

Then, when the wave goes through one cycle starting from the midpoint in 
the second medium, the beam transfer matrix becomes 



Mi 



e"^ 2/2 \ /cosh(r?/2) sinh(r//2) \ /e - ** 1 
e ^ 2/2 J Uinh(r//2) cosh(^/2)y V 



/ cosh(r?/2) -sinh(r]/2)\ /e~^ 2 / 2 On 
X 1 - sinh(r//2) cosh (77/2) J \ e ^/2 J ■ 



This arrangement of the matrices is illustrated in Fig. [2] 

The Mi matrix Eq. fl55|) contains complex numbers, but we are interested 
in carrying out calculations with real matrices. This can be done by means 
of a conjugate or similarity transformation [JO]. Let us next consider the 
matrix 

1/1 1\ /l i\_ 1 / e^/ 4 e^/ 4 



^"2 1-1 lJu J"7!Ue-^ 4 e— / 4 J- (56) 
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This matrix and its inverse can be written as 



We can then consider the conjugate transform of the M\ matrix 

M 2 = C Mi C-\ (58) 



with 

/cos(0 2 /2) -sin(0 2 /2)\ /e"/ 2 \ 
iKJ2 Un(0 2 /2) cos(0 2 /2) )\ e""/ 2 J 



' cos 4>i — sin 0i 
v sin 0i cos 0i 



3""/ 2 \ fcos(0 2 /2) -sin(0 2 /2)\ 
e^V Uin(0 2 /2) cos(0 2 /2) J' {0J> 

The conjugate transformation of Eq. fl58l changes the boundary matrix B(rj) 
of Eq. (|5ip to a squeeze matrix 

S fo)=( C r e -/ 2 )' (60) 

and the phase-shift matrices -P(0i) of Eq. (l52j) and Eq. (154"]) to rotation ma- 
trices 

R(24> i ) = ( C °fj} ~^), (61) 
v w V sm(0i) cos(0i J ' v ; 

with i = 1, 2. 
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